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5 ' Abstract 

g^ , 

We demonstrate that a simple model successfully describing experimental data 
Qh! for the process e~^e~ — > in can also qualitatively account for the data of CLEOII 

pi i' and ALEPH obtained recently for the decay r — > AnVr- The model is based on 

^ ' the assumption of the ai(1260)7r and cjtt dominance as intermediate states. Our 

observation is in contrast with the claim by ARGUS that the /o-meson signal in the 
r — > At^Vt decay can not be explained by the aivr intermediate state. 
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1 Introduction 



At present hadron processes at low energies (below J/ ip ) can not be satisfactorily described 
by QCD. Therefore, in order to understand the phenomena in this energy range it is 
necessary to apply different phenomenological models and compare their predictions with 
the accessible experimental data. Investigation of exclusive channels in e^e~ annihilation 
into hadrons at low energies and semihadronic decays of r meson is very attractive for this 
purpose. It provides the important information on the bound states of light quarks and 
elucidates the role of different mechanisms in these processes. 

Production of four pions is one of the dominant processes of e^e~ annihilation into 
hadrons in the energy range from 1.05 to 2.5 GeV. Due to the conservation of vector current 
(CVC) the cross section of this process is related to the probability of r — > AtiVt- decay 
Therefore, all realistic models describing the first process, should also be appropriate for 
the description of the other one. 

One of the main difficulties in the study of four pion production is caused by the 
existence of different intermediate states via which the final state could be produced, such 
as WK, pa, ai(1260)7r, /ii(1170)7r, p^p~, a2(1320)7r, 7r(1300)7r. The relative contributions 
of the mentioned processes can't be obtained without the detailed analysis of the process 
dynamics. Some information on the process e~^e~ Air has been obtained in [2-8], where 
the investigation of e~^e~ annihilation into hadrons was restricted by measurements of the 
cross sections only. However, the abundance of various possible mechanisms and their 
complicated interference results in the necessity of simultaneous analysis of two possible 
final states (27r+27r^ and 7r+7r^27r°) which requires a general purpose detector capable of 
measuring energies and angles of both charged and neutral particles. 

In the energy range below 1.4 GeV the detector of this type CMD-2 is operating at 
VEPP-2M collider in Novosibirsk 0. During a few last years very large data samples 
have been obtained which open qualitatively new possibilities for the investigation of the 
multihadronic production in e~^e~ annihilation. 



Very recently, a model-dependent analysis was performed ||T^ of both possible channels 
in e~^e~ annihilation into four pions at energies 1.05-1.38 GeV, using the data collected 
with the CMD-2 detector. The discussion of previous works devoted to this theme is also 
present in [|^. 

The detailed analysis of the process e^e~ — > Air unambiguously demonstrated that the 
main contribution to the cross section in the energy range 1.05 - 1.38 GeV , in addition 
to previously well-studied a;7r° [|TT], [T^, is given by the pirn intermediate state. More- 



over, the latter is completely saturated by the aivr mechanism. The contribution of other 
intermediate states was estimated to be less than 15 % . 

In this paper we use the assumption of the aivr dominance for comparison of the available 
experimental data for r —>■ AtcUj. decay |T^, |T^, |T3], with the prediction of our model 



successfully described the data of e~^e An. 
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2 Results and discussion 



The initial hadron state (referred to as p), which decays into four pions has the p-meson 
quantum numbers. Due to the conservation of vector current the probabihty dTi and 
of T~ — > 71~it^71~it^i't and r~ n'TT^Tc^TT^Ur decays, respectively, can be written as 

^II = %^(,ni + 2.)(,„? - (1) 

ds 967r3m3 ^ ^ ^ ^ W1 + W2 

where R^j^ is the ratio of the cross section e'^e~ Att and e+e^ l^^f^~, dWi and dW2 
are the probabilities of p~ decays into 7r^7r+7r^7r° and 7r^7r°7r°7r*^, respectively. Let dW^, 
and dWi be the probabilities of p° decays into -k^-k^tx^-k^ and Tr+Tr^vr'^vr^, then due to the 
isospin invariance, we have 

W^ = ]^W^ + W^ , W2 = \w^. (2) 

The explicit forms of the matrix elements, corresponding to PVj, are presented in the Ap- 
pendix. In order to get the predictions for r decay we neglect the interference between ujtt 
amplitude and aivr, and rewrite (|l|) in the following form: 



dVi cos' 9 



ds 967r^m^ 

~T~ = nr 3 3 + 2s - s) R27.+27.' -rrr- , (3) 

where dW^ is the probability of p~ — > cjtt" decay, R^,, is the ratio of the cross section 
e+e" ijjTi and e+e~ — > p^fi~- 

To get R2-K+2-K- and Ri^T, we used the cross section of e+e~ annihilation in the energy 
range 1.-2. GeV obtained by several experimental groups and represented in Figs. |I] and 
^ together with our fit. Our approximation of the experimental data took into account 
possible uncertainties in the data. Assuming that the total cross section e^e~ 27r+27r^ 
is saturated by the Oivr mechanism, we calculated the cross section e"'"e~ 7r+7r^27r°. The 
comparison of our approximation with the experimental data is shown in Fig. ^. One can see 
some difference between our expectation and data at energies above 1.6 GeV. We suppose 
that this difference can be explained by the possible systematic errors or by contributions of 
another intermediate states which become more important at higher energies. However, we 
expect that this difference does not affect very much our predictions for r since this energy 
region has the additional suppression factor due to limited phase space in case of r decay 
(see (|1|)). In order to check this we compared the distribution over four-pion invariant mass 
with the recent data obtained by CLEOII |ll6|] and obtained that our interpolation of e+e~ 



'ml + 2s) (m^ - sf 



T 



' dW^ , „ dW^ 

R..^ + R2.^2n-^ 



data was in acceptable agreement with the r decay data. 

In order to fix the parameters of our model we have used the data of e~^e~ — » 47r p] 



and r — >■ 2'7r°7r 2//^- [Q. Since the data in [|T^ were obtained below the threshold of 



OiTT creation it was difficult to extract independently both the mass and width of Oi. The 
mass of ai was taken from the PDG table and the width was obtained as a result of 
optimal description of e~^e^ 11^71^271^ data. Here we use the same ai mass and the 
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width we get as a result of optimal description by our matrix element (see Appendix) of 
three pion invariant mass distribution in r~ — >■ 27r°7r~^'T- decay |T^. We used this way since 



the data of r decay are more sensitive to the parameters of Oi meson. The comparison 
of our description of three pion invariant mass distribution with the data is presented 
in Fig. ^. We obtained the value of oi width which also provides a good description of 
e^e~ 4:71 data (see Figs. |], H). In [|l8| it was obtained the evidence that ai meson has 



significant probability to decay into three pions through air intermediate state. The data 
analysis of e~^e~ An also confirmed this statement. Therefore, here we take into account 
the admixture of (ttt to the ai decay amplitude and the parameters of this admixture we 
extracted from e~^e~ Att data. For comparison we also present the results obtained 
without arc contribution taken into account. 

When we get the parameters of our model, we can pass to the description of 
Ti^-K^-K^-K^Vj. decay. The most interesting information on the mechanism of four pion 
channel can be obtained from two-pion mass distributions. Fig. shows the comparison of 



our predictions with the data of CLEOII detector obtained without subtraction of ujn 



contribution. One can see a rather good agreement with this data. The same comparison 



with the ARGUS data [|I^ is made in Fig. |g. In this case our predictions are not in a 
rough contradiction with the data, but the agreement is essentially worse than with that 
of CLEOII, especially in the invariant mass range where the enhancement due to p meson 
can be seen. 

The data obtained with the subtraction of ujti~ contribution allows one to make a more 
detail comparison of the differential distributions predicted within the assumption on Oivr 
dominance. For this purpose we have used the data obtained by ALEPH (see Fig. 0) 



and very recent high-statistics data of CLEOII [|l6l (see Fig. |^). In the first case we see a 
good agreement although we can not take into account the detector efficiency and energy 
resolution. The data of CLEOII, which have essentially higher statistics, the agreement 
is a little bit worse. Unfortunately, in this data the contributions of background events 
(such as K^TT~i^^Vj. and K^tt^t[~tt^Vj. ) were not subtracted, though their fraction was 
significant (about 8%). The biggest discrepancy can be seen for the 7r"''7r~ invariant mass 
distribution. However the data of CLEOII in this case systematically differ from data 
obtained by ALEPH (see Fig. |c). 

In spite of some disagreement between the data of different groups we can conclude 
that the assumption of ai dominance is in qualitative agreement with all available data. 
The quantitative comparison of our model can be made only by full simulation which takes 
into account the energy resolution and detector efficiency as well as all possible background 
sources. On the base of our analysis we come to conclusion that the previously obtained 
upper limit (less than 44% at 95% CL) [|l3l for ainp^ contribution to the r AtiVt- decay 
is hardly correct. Note that in some theoretical works based on the effective chiral theory 
of mesons [|l^] the predictions of Oivr dominance were made but without detail comparison 
with experimental data. 
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3 Appendix: A model for e^e~ An and r AnUj 
processes 



To describe four pion production we use a simple model assuming quasitwoparticle 
intermediate states and taking into account the important effects of the identity of the final 
pions as well as the interference of all possible amplitudes. We use the notation for the 
polarization vector of the initial hadron state ( p), superscript corresponds to isospin. Here 
we consider the contributions of ai(1260)7r and ljtt) intermediate states to the amplitudes 
p ^ 47r in the energy range under study. Other contributions to these amplitudes from 
broad resonances having the masses close to the threshold of pir production (02 (1320), 



7r(1300) etc.) are discussed in |^ 



As known, the form of the propagators 1/ D{q) is very important for analysing the data. 



It was found in that it is necessary to take into account the dependence of the imaginary 
part of the propagators (width) on virtuality while the corresponding corrections to the 
real part which can be expressed through the imaginary part by dispersion relations are 
not so important. We represent the function D[q) in the form 

D(g)=g^-M^ + .Mr^, (4) 

where M and F are the mass and width of the corresponding particle, and the function 
g{s) describes the dependence of the width on the virtuality. If = then D = iMT, 
in accordance with the usual definition of mass and width of the resonance. In the case of 
the p-meson we used the following representation for the function gp{q): 

gp{s) = s-'/\s - 4mr^\ (5) 

where m is the pion mass. 

We include into consideration two the most important channels of ai Svr decay 
: ai ^ pn ^ Stt and ai ^ an ^ Sir. Taking into account the quantum numbers 
of pion and ai-meson, we can write the matrix elements corresponding to the processes 
p{P) -> ai{q)n{p), ai{q) p{P')7i{p), and ai(g) cr(P')7r(p) as 

Tip^a^n) = F-^a,^e'^'%P^el- P,~el)q,A':ct>'\ (6) 
T{a,^pn) = F,,,^e'^\,Al{P'/: - Ple';W\ 
Tia^^an) = Fa,.^{q,At - q,A';^)P'^p,r* , 

where a, b, c are isospin indices, A'^ and are the polarization vectors of Oi- and p-mesons, 
(p^ is the pion wave function, Fa^p^, , Fa^^n and Fp^^^ are form factors depending on the 
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virtuality of initial and final particles. The explicit form of these form factors in the energy 
range considered is not very essential (it contributes to the theoretical uncertainty of the 
model). The matrix element of the transition T(p — ^ tttt) and T{a tttt) reads 

T{a^nn) = F^^^r*^* , (7) 

where p^^^ and p^^^ are 4- momenta of the corresponding pions. The matrix element T(p ^ 
An) in the rest frame of p can be written as T = eJ , where e is the polarization vector of 
p. The corresponding probability of p decay is given by the usual formula 

,w(s)Jiii2.rsi'Hp.-P)n^,, (8) 

where P is the initial 4-momentum oi p { = E = ^/s , P = 0), = (ej,Pj) are 
the momenta of pions. Note that due to the helicity conservation in e^e~ annihilation 
only transverse space components (with respect to electron and positron momenta) of the 
4-vector are not zero. Using (^) and (|^ we obtain the following expressions for the 
contributions of the ai(1260)-meson to the current Ja^ = ia^^p-w + Joi-xttt in the process 
pO _^ 7r+(pi)7r+(p2)vr"(p3)7r-(j94) : 

K^-^fm = G[ti{pi,p2,P3,PA) +tl{pi,P4,P3,P2) +tl{p2,PuP3,P4) 
+tl(p2, P4, P3,Pl)+ tl(pi, P2, P4, P-i) + tl(pi, P3, Pi, P2) 
+ tl(p2,Pl,P4,P3) +tl(p2,P3,P4,Pl)] , (9) 

J^^a"; = ^[12(^1,^2,^3,^4) - t2(pl,P4,P3,P2) +t2(p2,Pl,P3,P4) 
-t2(p2,P4,P3,Pl) +t2(pi,P2,P4,P3) " t2 (Pl , P3, P4, P2) 
+ tl(p2,Pl,P4,P3) - t2(P2,P3,P4,Pl)] , 



where 



ti Pi,P2,P3,P4 = n V7T7 ^ \^ ^0 

P>aAP -Pi)P>p\Pi+P-i) 

x{(E-£4)[pi(E£3 -P4P3) -p^{Eei -PiPi)] 

-P4[^l(P4P3) -£3(^4^1)]}, 



t2(pi,P2,P3,P4) = n (TD \n r ^ V ^ 

Da^P - Pi)D^[pi + P3) 

x(P-p4)2[(E-e4)p2 + e2P4] 

Here \/DA{q) , \/Dij{q) and 1/D„{q) are propagators of ai , p and a mesons, Fa-^{q) is 
the form factor, G is some constant, z is the dimensionless complex constant. Similarly 
to (P), we obtain for the contributions of ai(1260) to the current Jai in the process p° — > 
7r+(pi)7r-(p4)7r°(p2)vr°(p3): 

KA^% = G'[-t„,(p4,P2,P3,Pl) +tai(Pl,P2,P3,P4) (H) 
- tai(p4,P3,P2,Pl) +tai(pi,P3,P2,P4)] , 
Jai^°^ = G 1*2(^2, Pi, P3,P4) - t2(P2,P4,P3,Pl) • 
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The contributions of the ai(1260)-meson to the currents J^^ in the process 
p~ 7r+(pi)7r°(p2)7r"(p3)7r"(p4) reads : 



j+O— 



-G'[ti(pi,P3,P4,P2) +ti(pi,P4,P3,P2) +ti(p2,Pl,P3,P4) 
+ tl(p2,Pl,P4,P3) - tl(pi,P3,P2,P4) - tl(pi,P4,P2,P3)] , 



K^^a^K = <^[t2(Pl,P3,P4,P2) +t2(Pl,P4,P3,P2)- 
-t2(Pl,P2,P3,P4) - t2(pi,P2,P4,P3)] , 



and in the process p — > 7r°(pi)7r°(p2)7r°(p3)7r (^4) it reads : 

J°7_^p^ = G [ti(p4,Pl,P2,P3) + tl(p4,Pl,P3,P2) + tl(p4,P2,Pl,P3) 

+ tl(p4,P2,P3,Pl) +tl(p4,P3,Pl,P2) +tl(p4,P3,P2,Pl)] , 

Ja7-^a7r = ^ [t2(Pl, P4, ^2, Ps) + t2(Pl, P4, Ps, P2) + t2(p2, ^4, P3, Pl)" 
-t2(pi,P2,P3,P4) - t2(pi,P3,P2,P4) " t2 (l>2 , P3, Pi , P4)] • 



(12) 



(13) 



(14) 



The function (s) in the propagator of ai has the form: 

^2Pl - £lP2 , £2P3 - ^3P2 



9a^ KS) 



KiQ) 



Dp{pi + P2) Dp{p2 + P3) -Da(pi + P3) 



+ 



+ 



I 1 2 
\z\^s 



Pi 



+ 



P2 



+ 



P3 



Da{P2 + P3) D^{pi + P3) D^(pi + P2) 



X 



X 



dPl dp2 dp3 (^^^^ (pi + P2 + P3 - <?) 
2£i2£22£3(27r)5 



(15) 



where = \/s and Q = 0. The first term in the braces corresponds to the decay 
ai — > 7r'^7r~7r°, and the second one to ai — > Stt*^. The function 5'(t(s) in the propagator of a 
is equal to 

g,{s) ^ {I - Am" / sf'\ (16) 

As a form factor, we used the function F{q) = (1 + ml^/h?) /(I + g^/A^) with A ~ 1 GeV. 
We found that in the energy region under discussion the amphtudes are not very sensitive 
to a value of A. 



The amplitude of the process p{V) — > uj{q)T^{p) has the form: 



;i7) 



where ep is the polarization vector of the tu-meson. The matrix element of the transition 
u ^ pTi can be written in the similar form. The c<j-meson contributes to channels p° — > 
7r+(pi)7r°(j92)7r°(p3)7r~(p4), 7r+(pi)7r"(p2)7r"(p3)7r°(p4) andp+ 7r"(pi)7r+(p2)7r+(p3)7r° 

The corresponding current is equal to 



Jo; = [ta;(P2,P4,Pl,P3) - ta,(p2,Pl,P4,P3) 

- ta;(p2,P3,Pl,P4)] + (P2 ^ P3) , 



(18) 
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where 



X{(£4P3 - £3P4)(PlP2) - P2(^4PiP3 " ^sPlPd 
-£2[P3(PlP4) - P4(PlP3)]} > 

£j is the energy of the corresponding pion, Fi^{q) is the form factor. Since the width of the 
cu is small, we set gui{s) = 1 in the propagator D^{q). 
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Figure 1: Energy dependence of the e+e — > 2tt^27t cross section 
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Figure 3: Energy dependence of the e^e tt+tt 2t[^ cross section 
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Figure 4: Three pion invariant mass distribution in the 7r~27r''i/T- decay. Blue curve 

is obtained taking into account ai — >• an decay. Red one corresponds only to ai pvr 
decay. 
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Figure 5: Distributions over invariant mass Mj„„(7r^7r°), Mrecoiz(7i"^), ^mt)(7r^7r^), 
Minvij^^T^^) for e^e^ -r^h^Ih^ process after a;7r° events subtraction. Blue curve is 
obtained taking into account a\ on decay. Red one corresponds only to a\ pit decay. 
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Figure 6: Distributions over invariant mass Mj„t,(7r+7r~), Mj„^(7r=''7r=''), Mrecoui'^'^) , and 
the distribution over the angle between the momenta of vr^ and vr" for e^e^ 27r+27r^ 
process. Blue curve is obtained taking into account ai an decay. Red one corresponds 
only to ai pn decay. 
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Figure 7: Distributions over invariant mass Mj, 



M(n"n) 



Minvin 7r°), and Mj„„(7r+7r ) 



for r — > 27r 7r"^7r°z/T- decay obtained by CLEO [|I^. Blue curve is the prediction obtained 
with tti an decay taken into account. Red one corresponds only to ai pn decay. 
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Figure 8: Distributions over invariant mass Mj„t,(7r+7r*'), Mj„^(7r~7r°), Mj„t,(7r+7r"), and 
Minvi'n'^'n'^) for 27i^7i^7T^Ur decay obtained by ARGUS |TB|. Blue curve is the 

prediction obtained with ai an decay taken into account. Red one corresponds only to 
Oi — >• pn decay. 
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Figure 9: Distributions over invariant mass Mj„t,(7r+7r'^), Mj„^(7r~7r°), Mj„t,(7r+7r"), and 
Minvi'n'^'n'^) for ^ 27r^7r+7r°z/^ decay after cjtt" events subtraction obtained by ALEPH 
\\14\\ . Blue curve is the prediction obtained with ai crvr decay taken into account. Red 
one corresponds only to ai — >• pn decay. 
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Figure 10: Distributions over invariant mass Mj„t,(7r+7r°), Mj„^(7r~7r°), and Mj„^(7r+7r^) for 
decay after ujt[~ events subtraction obtained by CLEO [T^ . Blue curve 
is the prediction obtained with ai an decay taken into account. Red one corresponds 
only to Oi — > pn decay. 
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